WMAP anomaly : Weak lensing in disguise 
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Statistical isotropy (SI) has been one of the simplifying assumptions in cosmological model build- 
ing. Experiments like WMAP and PLANCK are attempting to test this assumption by searching 
for specific signals in the Cosmic Microwave Background (CMB) two point correlation function. 
Modifications to this correlation function due to gravitational lensing by the large scale structure 
(LSS) surrounding us have been ignored in this context. Gravitational lensing will induce signals 
which mimic isotropy violation even in an isotropic universe. The signal detected in the Bipolar 
Spherical Harmonic (BipoSH) coefficients Afl^ by the WMAP team may be explained by accounting 
for the lensing modifications to these coefficients. Further the difference in the amplitude of the 
signal detected in the V-band and W-band maps can be explained by accounting for the differences 
in the designed angular sensitivity of the instrumental beams. The arguments presented in this 
article have crucial implications for SI violation studies. Constraining SI violation will only be pos- 
sible by complementing CMB data sets with all sky measurements of the large scale dark matter 
distribution. Till that time, the signal detected in the BipoSH coefficients from WMAP-7 could also 
be yet another suggested evidence of strong deviations from the standard ACDM cosmology based 
on homogeneous and isotropic FRW models. 
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The Cosmic Microwave Background (CMB) anisotropy 
measurements are one of the cleanest probes of cosmol- 
ogy. The CMB full sky temperature anisotropy mea- 
surements have been used to test the assumption of the 
isotropy of the universe. Ever since the release of first 
year data of the Wilkinson Microwave Anisotropy Probe 
(WMAP), statistical isotropy of the CMB anisotropy has 
attracted considerable attention. The study of full sky 
maps from the WMAP 5 year data [iHll and the very 
recent WMAP 7 year data [jl , has led to some intriguing 
anomalies that may be interpreted to indicate deviations 
from statistical isotropy. 

The CMB temperature anisotropics are assumed to 
be Gaussian, which is in good agreement with current 
CMB observations. Hence the two point correlation func- 
tion contains complete information about the underlying 
CMB temperature field. The two point correlation func- 
tion can be expressed in terms of the spherical harmonic 
coefficients of CMB temperature maps. 



function in the Lcgcndrc polynomial (P/) basis. 



C(fii,n2) = (Ar(ni)Ar(n2)) 
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In the case of statistical isotropy, the correlation function 
depends only on the angular separation between the two 
directions and not on the directions fii and 62 explicitly. 
This property makes it possible to expand the correlation 
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The angular power spectrum, C;, appearing in Eq. [3] 
encodes all information in the covariance matrix for the 
statistically isotropic case. This however is only true in 
the case of an unlensed CMB sky as will be evident from 
the primary message of the article. 

In the absence of statistical isotropy, the correlation 
function explicitly depends on the two directions fii and 
h-2. In this case, the covariance matrix has been shown 
to have non-vanishing off-diagonal elements. This feature 
of the covariance matrix is captured well in the Bipolar 
Spherical Harmonic (BipoSH) representation which was 
introduced by Hajian & Souradeep (HS) [1, @]- The Bi- 
poSH form a complete orthonormal basis for functions 
defined on 5^ x 5^. The CMB two point correlation 
function can be expanded in the BipoSH basis in the fol- 
lowing manner, 

C(ni,n2)- K''AYh{^^)®yiA^2)}LM (4) 
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LMI1I2 mim2 

where Cf"^ , ,„ are the Clebsch- Gordon coefficients, the 
indices of which satisfy the following relations, |/i — Z2I < 
L < li + I2 and toi -I- m2 = M . 

These BipoSH coefficients can be expressed in terms of 
the covariance matrix derived from observed CMB maps. 
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In the case of statistical isotropy the only non-vanishing 
BipoSH coefficients are A'^^ and they can be expressed in 
terms of the CMB angular power spectrum C; , 



(6) 



where IT; = \^2l + 1. We reiterate that the above dis- 
cussion is true only for the case of unlensed statistically 
isotropic CMB temperature fields. 

Weak lensing of the CMB photons due to scalar per- 
turbations leave measurable imprints on the CMB two 
point correlation function. The lensing modification to 
the CMB angular power spectra have been well studied 
0, Q ■ It is well known that lensing introduces coupling 
between different multipole moments of the temperature 
field, which otherwise are not expected to be present in 
an unlensed statistically isotropic CMB sky (See Eq. [3]). 
This coupling between the various multipole moments 
arises due to the coupling between the projected lensing 
potential ^(n) Q and the unlensed CMB temperature 
field. This feature again can be most generally captured 
in the BipoSH representation. 

In the discussion that follows we have assumed that the 
unlensed CMB temperature field is statistically isotropic 
and that the projected lensing potential and the CMB 
temperature field are uncorrelated. Lensing remaps the 
temperature field. To leading order terms in the deflec- 
tion field A, the lensed temperature field T(n) can be 
expressed in terms of the unlensed temperature field T, 



f(n) = r(n + A) 



(7) 



where the lensing deflection fleld on the sky can be ex- 
pressed in terms of the gradient of the projected lensing 
potential fleld on the sphere as, 



Aa = Va*(n). 



(8) 



The two point correlation of the lensed temperature field 
can be expressed as follows, 

(f(fii)f(n2)) = (T(ni + Ai)r(n2 + A2)) 

= (T(ni)r(n2)) (9) 

+ (V'^V(fii)VaT(fii)T(n2)) 

+ (V'^7A(ri2)Var(n2)T(ni)) + 0{i'^) . 

The corrections to the CMB angular power spectrum due 
to lensing arise only due to terms which are of 0{ip'^), 
whereas the corrections to the BipoSH coefficients can 
be shown to be only due to terms linear in the lensing 
deflection field 

Without making assumptions about the isotropy of the 
lensed temperature field, the two point correlation can 
be most generally expressed in terms of the BipoSH co- 
efficients and the spherical harmonic coefficients of the 
lensing deflection fleld ipim, 
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Since we have already assumed the unlensed temperature 
field to be isotropic, the above equation can be further 
simplified and expressed (Using Eq. |6] ) in terms of the 
unlensed CMB angular power spectrum, 

+ 2afJ/(^;™,Qi,Q.).(ll) 

Evaluating Eq. ^ in the harmonic space allows us to 
obtain the following expression for the coefficient CtfJ)^, 
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Ci,F{h,L,h) + Ci,F{h,LM) 
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where, 



F{h,LM) = 



[l2il2 + l) + L{L + l)-h{h + l)] 



Note that lensing by scalar density perturbations nec- 
essarily generates only the even parity (i.e. li + I2 + L 
is even) BipoSH coefflcients which is apparent from Eq. 
[T2] due to the presence of C;^g;^g which vanishes for odd 
parity {li +I2 + L \s odd). 

The above discussion of BipoSH coefficients is in terms 
of the HS estimator which differs from the estimator used 
by the WMAP team. The WMAP estimator is vahd only 
for even parity Bipolar coefflcients unlike the HS estima- 
tor. For even parity coefflcients, the two estimators are 
related by the following expression. 
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In discussions that follow we use the WMAP estimator 
for the BipoSH coefflcients. 

Motivated by the WMAP detections of isotropy viola- 
tion Q in the V-band and W-band maps in ecliptic coor- 
dinates, we study the corresponding BipoSH coefflcients 
Afl^ and ^^;_|_2 that arise due to lensing. These coeffl- 
cients take up an extremely simple form proportional to 
the angular power spectrum given by. 
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(14b) 



where corrects the CMB angular power spectrum C; 
for convolution of the instrumental response beam func- 
tion. 

We repeat the WMAP analysis to measure the Bi- 
poSH coefficients from the measured V-band and W- 
band maps. In our analysis we account for the effects of 
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FIG. 1: The amplitude of the signal detected in the 
BipoSH coefBcient Afj^ in the V-band map is less than 
the amplitude of the signal in the W-band map exactly 
corresponding to the expected difference due to the two 
instrumental beam widths. It is remarkable and 
intriguing that it is possible to explain the V-band and 
W-band detections with a common, consistent value for 
the quadrupolar component of the projected lensing 
potential 1^20 ■ See Table U for the beam specifications of 
WMAP and the best fit value of the parameter ip2o- 
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FIG. 2; The signal in the BipoSH coefficients Af°i_^_^ 
predicted due to lensing (See Eq. I14bl ) plotted against 
the detected signal from WMAP 7 year maps. The 
value of the parameter V'20 is set to the best fit value 
derived from the fit obtained to the BipoSH coefficients 
Af^ and can be found in Table U . 



beams. Since the harmonic coefficients are coordinate 
dependent quantities, we specifically mention that this 
value of tp20 inferred is in ecliptic coordinates. The re- 



tire anisotropic noise, instrument beam, foregrounds and 
masking, each of which can mimic signals of isotropy vi- 
olation. We estimate the bias introduced due to all the 
above mentioned effects (except for the non-circularity of 
the instrument beam) through simulations. This bias is 
then subtracted from the signal inferred from the analysis 
of the WMAP V-band and W-band maps to isolate the 
systematic effects. We do not perform the inverse noise 
weighting that results in a much larger estimate for the 
errors here and lowers the significance of the detections 
in our analysis. The broad features of the detections are 
however retained and suffice for the arguments presented 
in this article. More elaborate analysis incorporating in- 
verse noise weighting on harmonic coefficients aim of the 
temperature field and the non-circularity of the beam is 
underway and will be reported in a more detailed publi- 
cation. 

In our analysis we use the best-fit ACDM CMB an- 
gular power spectrum Cj generated using the publicly 
available Boltzmann code CAMB [l2|. We use the lens- 
ing field harmonic coefficient -020 which appears in Eq. 
I14al as a free parameter, as currently no measurements 
of the projected lensing potential exist at the largest an- 
gular scales. We perform a simple minimization tech- 
nique to estimate the bestfit value of this free parameter. 
Note that we carry out the fitting procedure only with 
the BipoSH coefficients Af^ , since the detections in these 
coefficients are highly significant. Our analysis yields a 
value of ip20 = 2.25 x 10~^ with which we are able to 
fit the detected signal in the BipoSH coefficients Afj^ in 
both V-band and W-band maps equally well, after ac- 
counting for the effects of the respective instrumental 





FWHM (in degrees) 


per d.o.f 


V'ao 


V-band 


0.326 


0.818 


2.25 X 10"^ 


W-band 


0.202 


0.854 



TABLE I: WMAP beam specifications from lLAMBPXl 
site. The per degree of freedom (d.o.f) for the 
best-fit "020 parameter that is obtained by fitting Eq. 
[Tia] to the BipoSH coefficient Af^ detections in WMAP 
7 year maps. 

suits of the whole analysis are summarized in Table |T] and 
Figure [T] The exceedingly good fit to the data as sug- 
gested by the reduced values in the above table are due 
to over estimates of the standard deviation on the data 
points in absence of inverse noise weighting. The best fit 
value of the parameter il'20 is then used to predict the sig- 
nal due to lensing in the BipoSH coefficients Af'j_^_2 (See 
Eq. I14bl ) . The results are plotted against the detections 
found from WMAP 7 year maps and are depicted in Fig- 
ure [2] In this case the predicted signal does not seem to 
match the detections, however this could be because of 
other systematic effects which remain unaccounted, like 
the non-circular beam. Nevertheless it is interesting that 
the predicted signal for these BipoSH coefficients follows 
the trend seen in the detections particularly at low mul- 
tipoles. 

We have clearly established that gravitational lensing 
can introduce significant corrections to the BipoSH co- 
efficients rendering them non-zero even in an isotropic 
universe. We have also demonstrated that the difference 
in amplitude of the signal detected from V-band and W- 
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band maps can be explained by accounting for the re- 
spective beam angular sensitivities. 

The best fit ACDM cosmology predicts the quadrupo- 
lar power in the projected lensing potential power spec- 
trum Cj"'' to be of the order of ^ 10~^ requiring the 
corresponding harmonic coefficients 'ip2m to be of order 
~ 10"''. This suggest that the value of the quadrupole 
estimated in our analysis would be highly improbable. 
Hence in the realm of standard ACDM cosmology the ob- 
served signal in the BipoSH coefficients may not be com- 
pletely explained by accounting for modifications due to 
lensing. Hanson et. al. [3| have argued that the detec- 
tions found in the BipoSH coefficients may be explained 
by accounting for the non-circular nature of the WMAP 
instrumental beam. The detections of the BipoSH co- 
efficients in the WMAP 7 year data may be completely 
explained by accommodating the total contribution from 
both these effects. 

After correcting for the systematic effects of the beam, 
it will be important to account for the lensing modifica- 
tions to the CMB two point correlation function as has 
been discussed throughout this article. It is also perti- 
nent to note that odd parity BipoSH of HS that may be 
associated with the effect non-circular beam with a spe- 
cific scan strategy could possibly distinguish between the 
two effects. Otherwise, this would require a completely 
independent measurement of the lensing potential field 
and can be only provided by LSS measurements. In the 
following section we discuss the LSS measurements that 
will be essential to estimate the projected lensing poten- 
tial 

The projected lensing potential is constructed by form- 
ing a linear weighted sum of the gravitational potentials 
along the line of sight. Since the dominant contribution 
to the gravitational potentials is due to dark matter, an 
all sky map of the large scale distribution of dark mat- 
ter will be needed. This will be made possible through 
gravitational lensing surveys. A significant contribution 
to the projected lensing potential power spectra comes 
from LSS below redshift z < 5 [0]. However to esti- 
mate the low multipole moments of the projected lensing 
potential, it will suffice to map the dark matter distri- 
bution upto redshift z < 2. This will allow a reasonable 
estimate of the lensing contribution to the detections in 
the BipoSH coefhcients. 

The Cosmic Evolution Survev [UOSMOSI has mapped 
the dark matter distribution out upto redshift z — 5, 
however only on a small patch of the sky covering 2 
square degrees. Similar deep surveys with much larger 
sky coverage will be required to achieve this goal. This 
will be made partially possible with upcoming surveys 
like Large Synoptic Survey Telescope LSS T, Dark En- 



ergy Survey [DE5] and lUnLIDl 

Next, we discuss the possibility of completely explain- 
ing the detection of the BipoSH coefficients, as arising 
from the lensing modifications alone. The statistically 
significant detections of the quadrupolar bipolar power 
spectrum in the recent WMAP 7 year results [1] can be 
explained by accounting for the lensing corrections to the 
BipoSH coefficients. This explanation would imply that 
the large scale distribution of dark matter surrounding us 
happens to have an extremely high quadrupole moment. 

There have been other anomalies, like the cold spot 
observed in the CMB maps, which may not be explained 
in the realm of the standard ACDM cosmology. It has 
been suggested that this anomaly in the CMB maps may 
be explained by correcting the observed CMB maps for 
integrated Sachs- Wolfe (ISW) effect due to an immensely 
large void in the LSS surrounding us [l^.[l5j. Further it 
has been argued that the presence of such large voids may 
only be explained by invoking non-Gaussianity in the pri- 
mordial density fluctuations [l^. The best fit value of 
'ip2a obtained in our analysis may not seem as improba- 
ble if specific non-Gaussian initial conditions are invoked. 
More follow up work is motivated by our results. 

We have demonstrated that the difference in the am- 
plitude of the signal detected in WMAP V-band and 
W-band maps can be explained by accounting for the 
beam angular sensitivities. Further, Grocneboom et. al. 
p7| assess that the non-circular beam cannot explain the 
non-zero detections of the BipoSH coefficients, contest- 
ing the claims made by Hanson et. al. These argu- 
ments seem to allow for the possibility of the detection 
in the BipoSH coefficients having a cosmological origin. 
A weak violation of isotropy may result in a relatively 
large value for the quadrupole of the projected lensing 
potential which in effect could magnify the SI violating 
signal through lensing. In such a case the detections may 
actually suggest a violation of statistical isotropy |18l - [2l| . 

Finally, the value of the parameter ■i/'2o inferred could 
be interpreted as first explicit measurement of projected 
lensing potential at the largest angular scales. Gravi- 
tational lensing has been proposed as one of the tests 
for modified gravity models [22[, suggesting yet another 
possibility of constraining these models through measure- 
ments of the lensed BipoSH coefficients. 
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